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Abstract 

CIh' The first part of this paper contains an introduction to Bell inequalities and Tsirelson's 

theorem for the non-specialist. The next part gives an explicit optimum construction 
C ' for the "hard" part of Tsirelson's theorem. In the final part we describe how upper 

bounds on the maximal quantum violation of Bell inequalities can be obtained by an 
extension of Tsirelson's theorem, and survey very recent results on how exact bounds 
may be obtained by solving an infinite series of semidefinite programs. 



Keywords: Quantum information. Quantum computation. Quantum correlation. Bell 
' inequality, Non-local games 

00 ■ 

00 

: 1 Introduction 



Correlations between observed measurements are a fundamental resource in quantum in- 
formation. In his seminal 1964 paper. Bell demonstrated an inequality that must be 
satisfied by correlations obtained classically, but may be violated by a quantum correlation 
experiment. Subsequently, the term Bell inequalities has come to describe the set of inequal- 
^ ' ities that characterize correlations that can be obtained between events that can be well 

■ described by classical physics. Since the 1980s, a series of experimental results have been 

published that apparently demonstrate that quantum correlations do in fact violate Bell in- 
equalities |T2] . This raises the question of whether a "good" mathematical characterization 
of quantum correlation vectors can be obtained. 

The first part of the paper is intended for the non-specialist. We begin by giving a sur- 
vey of classical correlations, showing that they can be characterized using the well studied 
subject of L^-embeddings. Next we consider quantum correlations and give a statement 
of Tsirelson's theorem [5], which gives a partial characterization for the 2-party case. In- 
terestingly, as we show, this characterization can be restated in terms of L^-embeddings. 
Tsirelson's theorem has found many applications in quantum information, but a direct con- 
structive proof of his theorem is not readily available. Although this theorem is powerful, 
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even in the two-party case it does not solve the characterization problem completely. We 
discuss this, and give an additional necessary condition based on no-signalling. In many 
applications of Tsirelson's theorem it is necessary to construct a quantum realization for a 
given correlation vector. An explicit construction is not given in his paper, but we give one 
here. We begin by giving a simple construction for the two dimensional case of Tsirelson's 
theorem, which uses a two qubit state. This includes introducing many of the ideas needed 
for the general construction. 

In the second part of the paper, we give a new general construction, that is, states and 
operators which can represent all correlation functions. Our method is related to the theory 
of stabilizer states. We also demonstrate the optimality of this result, and give a scheme 
for a low dimensional approximation of an optimal representation. 

In the third part of the paper we discuss the problem of finding the maximum violation 
of a Bell inequality. We see that semidefinite programming (SDP) plays a crucial role in 
this optimization problem, and relate it to the results given earlier in the paper. Finally 
we review some very recent work involving an infinite hierarchy of SDPs that gives a non- 
polynomial time method of characterizing Bell inequalities and finding maximum violations. 

2 Classical correlations 

Let Ai, ...,An be a collection of n 0/1 valued random variables that belong to a common 
joint probability distribution. For 1 < i < j < n, we define new random variables AiAAj 
that are one when Ai = Aj and zero otherwise. Denote by {A) the expected value of a 
random variable A. The full correlation vector x based on ^i, An is the vector of length 
N = n + n{n — 1) /2 given by the expected values: 

X = {{Ai) , {A,AA,)) ^ ((^.)i<,<„ , (A,AA,)i<^<^.< J. 

Note that each element of the above vector lies between zero and one. Now consider any 
vector X = (xi, a;„, a;i2, G [0,1]^ indexed as above, which we will call an 
outcome. We consider two related computational questions: 

Recognition. When is an outcome x a full correlation vector? 

Optimization. For any c G what is the maximum value of (Fx over all possible full 
correlation vectors x? 

It turns out that the recognition problem is NP-complete, and the optimization problem 
is NP-hard. This follows from the fact that the set of full correlation vectors is in fact the 
cut polytope CUTn+i defined on the complete graph Kn+i- This polytope is defined as 
the convex hull of the 2^ full correlation vectors obtained by deterministically setting each 
random variable Ai to either zero or one. For details of the above and other facts about cut 
polytopes, see the book by Deza and Laurent [10]. For a vector u = (ui, Ud) the Li-norm 
of u is given by ||ti||i = X]f=i We have the following well-known characterization of the 
cut polytope. 

Li-characterization of full correlation vectors. 

The following two statements are equivalent: 

• An outcome x G [0, 1]^ is a full correlation vector. 



• There exist vectors G R'^, \ <i<n, d<N, for which 

Xi = IK^IIi, Xij = ||n* — n-'lli, ^ < i < 3 < n. 

Full correlation vectors provide an adequate model for correlations obtained in physical 
experiments at the classical level. Let us call the random variables observables. For example, 
with n = 3, Ai,A2,A3 could obtain the value one if a given McGill student has blond hair, 
weighs more than 80 kg or is more than 180cm high, respectively. We could obtain a full 
correlation vector by determining these three observables for all McGill students. 

In a quantum setting, things are very different. Firstly, it is difficult to apply the above 
model directly since at the quantum level it may be impossible to measure directly different 
observables for a given particle. Therefore the above model is replaced by a bipartite 
setting where the 0/1 random variables (observables) are labelled Ai, ...,Am and Bi, ...,Bn 
respectively. 

The (bipartite) correlation vector x based on random variables Ai, Am and Bi, Bn 
is the vector of length M = m + n + mn given by the expected values: 

X = {{Ai),{Bj),{AiABj)) 

= ((^«)l<i<m ' (-^i)l<j<n ' (^«^-^j)l<i<m,l<j<n)- 

As we will be concerned only with the bipartite case, we will simply use the term 
correlation vector where no confusion arises. As before, we call any vector x G [0, 1]*'''^ 
indexed as in ([T|) an outcome. Again we may define a polytope by considering the convex 
hull of the 2™''^" correlation vectors formed by letting each of the m + n random variables 
take value either zero or one. This polytope is called the Bell polytope Bfn,n and was 
apparently first considered by Froissart [13]. It turns out the membership and optimization 
problems given above are still NP-complete and NP-hard respectively (for references, see, 
e.g., [2] ). The characterization theorem generalizes in a natural way. 

Li-characterization of bipartite correlation vectors. 

The following two statements are equivalent: 

• An outcome x G [0, 1]*''^ is a bipartite correlation vector. 

• There exist vectors u^,v^ G W^, 1 < i < m, 1 < j < n, d < M , for which 

= Xm+j = ll^^-'lll, Xij = - U-'lli. 

The Bell polytope has been much studied. Valid inequalities for the Bm,n are often called 
Bell inequalities, although here we will reserve this term for the facets of Bm^n- These 
inequalities have been studied by many researchers, see for example [8], [22], [29]. The 
CHSH inequality is the only non-trivial facet of B{2,2) and is given by 

(^lA^i) - {A1AB2) - {A2AB1) - {A2AB2) < 

or equivalently 

xu - X12 - X21 - X22 < 0. (2) 

Although few Bell inequalities were known until recently, much is known about facets of 
the cut polytope, including several large classes of facets. In [1] a method is given to 



generate Bell inequalities from facets of the cut polytope, producing a large number of new 
inequivalent Bell inequalities. 

The correlation vector with m = 2, n = 2 given by 



clearly violates the CHSH inequality ([2]), so it follows there is no joint distribution function 
for the four random variables. This correlation vector cannot arise as the result of an exper- 
iment for which the rules of classical physics apply. An outstanding prediction of quantum 
theory, apparently confirmed by numerous experiments, is that this correlation vector can 
arise from observations at the quantum level. This fact has lead to many surprising appli- 
cations in quantum information theory, see for example Cleve et al.f^. It raises the issue 
of whether there is a good characterization of such quantum correlation vectors, the topic 
of the rest of the paper. 

3 Quantum correlations 

The postulates of quantum theory give a complete statistical description of the outcome 
of experiments at the quantum level. A two party quantum correlation experiment can 
be described by a quantum state and set of observables Ai, A^, Bi, ...Bn on a bipartite 
Hilbert space. It is assumed the two parties are spatially separated and that the observations 
are performed essentially simultaneously, so that there is insufficient time for the parties 
to communicate. For a given experimental outcome, the vector x defined by ([T]) is called 
a quantum correlation vector. The description given by the postulates does not appear to 
provide any tractable method to answer the recognition, optimization and characterization 
questions when applied to quantum correlation vectors. Such answers are provided, however, 
for one important case by a theorem of Tsirelson. A quantum correlation function is a vector 
y G ii"^" defined by taking the last mn coordinates of a quantum correlation vector, i.e.. 



1 1 1 1 2 + ^2 2-^/2 2-\/2 2-\/2 
2' 2' 2' 2' 4 ' 4 ' 4 ' 4 



(3) 




(4) 



Tsirelson's Theorem (0/1 version) [5j [27J. 

The following three statements are equivalent: 



. y = {{A,ABj)) e [0,1] 



mn 



is a quantum correlation function. 



• X = (1/2, 1/2, 1/2, (AiABj)) G [0, 1] is a quantum correlation vector. 

• There exist vectors u*, G R'^, 1 < i < m, I < j < n, d < m + n, for which 

. 1 . 1 . 

•^i — 11''^ II — X) ^m+j — 11^ II — '^ij — ll''^ II ■ 



where \\u\\ = u u. 



We call an experimental outcome unbiased if for all i and j we have {Ai) = (Bj) = 1/2, oth- 
erwise it is biased. A remarkable result implied by this theorem is that the recognition and 
optimization problems for correlation functions and unbiased quantum correlation vectors 



can be solved in polynomial time by semidefinite programming(SDP). Using the theorem, 
we can verify that ([3]) is a quantum correlation vector by exhibiting the vectors: 



1 _ I' -1 1 2 _ / 1 ^ ^^ 

" ^271' 271' 2^' "^271' 271' 

Furthermore, it can be verified by SDP that this is the maximum violation of ([2]), although in 
this case Tsirelson [5] has provided an analytic proof. These maximum quantum violations 
have many interesting applications, see e.g. |7|. The maximum quantum violation of any 
Bell inequality (like CHSH) that does not have terms involving the expectations {Ai) or (Bj) 
can likewise be found by using SDP. Unfortunately, most of the Bell inequalities produced 
recently [1] do not satisfy these conditions. For these inequalities the maximum quantum 
violation may only be achieved by a biased quantum correlation vector, and the above 
method cannot be directly applied. 

Tsirelson's theorem may not hold for experimental outcomes that are biased. Consider 
the outcome for m = n = 1 given by x = (3/4,3/4,3/4). If we set = (\/3/4, 3/4) and 

= {-^/3/4,3/A) then 

Xl = ||n*||, X2 = \\v^\\, Xl2 = — v^W, 

and the corresponding vector y = (3/4) is obviously a quantum correlation function. How- 
ever X is not a quantum correlation vector because it violates the no-signalling condition. 
This condition derives from the fact that the expectations {Ai) ,1 < i < m should be the 
same regardless of which measurement j the other party decides to make, due to the spatial 
separation of the two parties. Similar conditions should hold for the expectations {Bj). It 
is shown in [2] that a vector x satisfies the no-signalling condition if and only if it belongs 
to the rooted semimetric polytope defined by the inequalities: 

Xi ~\~ Xj ~\~ Xij ^2, Xi ~\~ Xj ^ij — ^? 
Xi j I '^ij ^5 ~^ j I '^ij ^* 

It is easy to see that unbiased quantum correlation vectors satisfy the no-signalling condi- 
tion. However, the vector x = (3/4,3/4,3/4) violates the first of these inequalities. 

It is tempting to conjecture that an outcome x is a quantum correlation vector if it satis- 
fies the no-signalling conditions ([6]) and the corresponding vector y is a quantum correlation 
function. However, consider the vectors 

_1 1 112 + 722-^/22-^/22-^/2 
" ^171' 171' 2' 2' 

_ 2 + ^/2 2-V2 

The outcome x satisfies ([U]), and y is a quantum correlation function, as shown by the 
vectors given in ([5|). Nevertheless, it is proved in [3] that x is not a quantum correlation 
vector. Perhaps even more surprising is a non-quantum outcome they exhibit for the case 
m = n = 3: 



1 2 

-, 1 < « < 6, xii = X22 = 0, Xij = - 

for all other 1 < i < j < 3. 



This gives an outcome x which satisfies ^ and for which the corresponding correlation 
function can even be obtained classically. For example with vectors 



= v^ = (0, 0, 0, 1/3), = v'^ = (0, 0, 1/3, 0), 
= (1/3,0,0,0), = (0,1/3,0,0) 

we have 

Xij = \\u^ — v^Wi 1 < « < J < 3. 

and can use the Li characterization theorem given in the previous section. 

Even though Tsirelson's theorem does not give a characterization of quantum correlation 
vectors, it can be extended to give a necessary condition that can be combined with the 
no-signalling condition. 

Necessary conditions for quantum correlation vectors f2]. 

li X = {{Ail , {Bj) , {AiABj)) G [0, 1]^^ is a quantum correlation vector then 

• X must satisfy the no-signalling conditions (0), and 

• There exist vectors u^, G R'^ , 1 < i < m, 1 < j < n, d < m + n, for which 

— 11^ II? ^m-irj — 1 1 1 1 5 '^ij — 11^ "^"^ 1 1 • 

where \\u\\ = u^u. 

Using this theorem, it can be shown that in the two previous examples the outcomes are not 
quantum correlation vectors. It also provides an efficient means of bounding the maximum 
quantum violation of general Bell inequalities by semidefinite programming (SDP). The 
strength of the no-signalling conditions with respect to various known Bell inequalities is 
discussed further inlCl 

It has recently been shown that the above conditions are not sufficient: in Doherty et 
al. [11] and Navascues et al. [21] a hierarchy of SDPs is presented which in some cases 
gives much tighter bounds. A characterization of quantum correlation vectors is given by 
an infinite hierarchy of conditions, each of which can be tested by solving an SDP. This 
work is reviewed in Section [51 It is still, however, an open problem to determine whether or 
not there is a polynomial time algorithm to determine whether a given vector is a quantum 
correlation vector. 

4 A proof of Tsirelson's theorem 

In this section we give an elementary description of the easy direction in Tsirelson's proof. 
For the proof, it is convenient to let the observables Ai, A^, Bi, Bn take values ±1 
rather than 0/1, and to consider the products {AiBj) rather than the differences {AiABj). 
An outcome is now given by 

X = {{Ai),{Bj),{AiB,)) 

l<i<m ' l<j<n ' *-^i)l<i<m,l<i<n)' 

and is called a quantum correlation vector if it can result from a quantum experiment. Simi- 
larly we redefine a quantum correlation function. In this section we use the ket-bra notation 



where \u) denotes a (possibly complex) vector, {u\ denotes the transpose of its complex con- 
jugate, and {u\v) denotes inner product. Using these new notations, the theorem takes the 
following equivalent form [5l [27] . 

Theorem 1 (Tsirelson's Theorem (±1 version)) For real m x n matrix Cij, the fal- 
lowings are equivalent: 

(1) y = {cij) £ [—1,1]"^^ is a quantum correlation function. 

(2) X = (0, 0, 0, Cjj ) G \^_\^\\rn+n+mn ^ quantum correlation vector. 

(3) There exist two sets of unit vectors {n*}^;^, {v^}^^i E W^, d < m + n, for which 

Cij = (n*|i;-') . 

(4) There exist two sets of vectors {u^^'^i and {v^Y-^^ on i^m^C™-'") such that they satisfy 

< 1, Iw-'l < 1 for all i and j, 

Cij = {u'\v^). (7) 

(5) There exist a state on a composite Hilbert space TLa ®'Hb and two sets of Hermitian 
operator {Ai}^^ on a Hilbert space Ha and {Bj}^^^ on a Hilbert space Hb such that 
\Ai\ < I, \Bj\ < I, and 

Cij = Tv{Ai®Bjp). (8) 

Proof: If (1) holds, some quantum correlation vector x' must be consistent with y. 
By switching the outcomes +1 and -1 we see that the vector x" formed by setting x'l = 
— x'j, i = 1, ...,m + n and otherwise setting x'l^ = x[ - is also a quantum correlation vector. 
The implication (1) ^ (2) then follows from the convexity of the set of quantum correlation 
vectors by setting x = {x' + x")/2. 

The implication (2) ^ (3) follows from the postulates of quantum theory. Indeed, cor- 
responding to the given quantum correlation vector there must exist observables Ai, Am 
on a Hilbert space TLa-, observables ...,i?n on a Hilbert space TLb-, and a pure quantum 
state given as a unit vector on TLa 'S'TLb, where denotes tensor (Kronecker) product. 
For 1 < i < m and I < j < n, let |a*) = Ai(g) Ib \tp) and \ V) = Ia® Bj Then |a*) and 
\V) are (possibly complex) unit vectors of length, say, t, such that (a*|6^) = (AiBj). We 
may replace them with real vectors and of length 2t by writing the real and complex 
coefficients as separate coordinates, maintaining the same values of the inner products. The 
set of m + n real vectors u^, u™, v^, have all the properties of part (2) except pos- 
sibly their dimension 2t > m + n. However the unit vectors lie in a subspace of dimension 
d < m + n, which preserves their inner products. We can also prove (5) ^ (4) in the same 
way. 

Now, we prove (4) =^ (5) and (3) =^ (1). Suppose the condition (4) (or the condition 
(3)) holds, and ^ = min ^dimspanju*}™^, dimspanjti-'}"^ (or ^ = d for the condition (3)). 

We choose a Hilbert space Ti whose dimension is 2'^ , where = I even ^, and = 
for odd ^. Then, there exists an irreducible representation of the Clifford algebra on TC, 
and we can choose a set of ^ Hermitian operators {^fc}|=i as an irreducible representation 
of a generator set of the Clifford algebra (see [A]); that is, satisfies the following 

anti-commutation relation: 

X,X, + XjXi = 26ijl. (9) 



Then, from the above relation, Hermitian operators {X^C^^Xk}^^^ commute with each other. 
Thus, they have simultaneous eigenvectors, and moreover, their corresponding eigenvalues 
are ±1. Therefore, there exists a state \^) € TCa 'S>'Hb which is a simultaneous eigenvector 
of {Xk satisfying 

Xk®Xkm = {-ir\^) (10) 

where is or 1. Then, and satisfies 



^{(-ir {^\XkXi^i\^) 



{^\Xk^Xi\^) 

+{-ir^ {^ixiXk^im 
= \{-ir{'^\{XkXi+xiXk)(^i\^) 
= i-irskii, (11) 

where we use Eq. (llOp in the first equality and Eq. ([9]) in the third equality. Finally, we define 
Ai and Bj as 

Ai = ^{u%Xk 

k=l 

B, = Y.{-ir{v")kXk, (12) 

k=l 

where u'* and v'^ are vectors derived by projecting ti* and onto span{n*}^j^ in the case of 
dimspanjti*}™^^ < dimspanju-^j^^j^ and onto span{f-'}™^j^ otherwise (for condition (3), we 
can just choose u'^ = and v'* = v*). By their definition, we can easily see that = |n'*p/ 
and -B? = \v'^\^I. This fact guarantees that \Ai\ < I and \Bj\ < I (for the condition (3), Ai 
and Bj take eigenvalues ±1). Finally, we can calculate {^\ A^ Bi\^) as follows: 

{^\A,^Bj\^) 

= Y.{-ir'{u%{v'^)i{^\Xk0Xi\^) 

kl 

= Y.^-lf''^iu'')k{v'')iSki 

kl 

= {u"\v'^) 

= {u'\v^) (13) 

where we use Eq. fllip in the second equality. □ 
In the above proof of (4) (5) and (3) ^ (1), in order to satisfy the equation 
(^'1 (E> -Bj 1^*) = (^u^\v^'j, Tsirelson used the simultaneous eigenvector 1^*) correspond- 
ing to eigenvalue "1" of {Xk}l^^ [26]. With Ai and Bj defined by Eq.((l2]) and a'' = 0, 
we showed that we can choose an arbitrary simultaneous eigenvector of {Xk ^ Xk}k=i by 
slightly changing the definition of Ai and Bj from Tsirelson's original definition. As we will 
see later, this modification of the proof is actually important for finding an explicit optimal 
construction of Ai and Bj . 



5 Representing two dimensional vectors 



In order to illustrate the concepts, we give a construction of a quantum setting that can 
realize the inner products of two dimensional vectors. Let u^, vP^ and f^, be m + n 
unit vectors in and let z = (21,^2)- Define 



Now define observables 



C{z) 



Zl Z2 
Z2 -Zl 



l,...,m j 



Let I*) = [l/\/2, 0, 0, l/VW which corresponds to the state '""^^^^^ . Then it is easy to 
verify that Eq. (jl3p holds. For example, with 

= [1,0]^, u^ = %lf, 

= [I/V2, 1/V2f, = [I/V2, -1/V2f 

we may verify, for instance, that 

Incidentally, this gives a quantum setting for the maximum violation of the CHSH inequality 
described in Section [21 when the inequality is restated in ±1 terms. 
Here is how the quantum setting is obtained. We start with 



1 

L -1 



1 

1 



It is easy to verify that the anti-commutation relation ([9]) holds (that is, Xi,X2 are a 
representation of the generators of the Clifford algebra C/(M^)) and that C{z) = X{z) as 
given by with = 0. To get the state 1"$) we first construct the operator A given by 



2A = Xi^Xi + X20 X2 



10 1 
0-110 
1-10 
10 1 



Using Maple we find that A has maximum eigenvalue 1 with corresponding eigenvector 
[1, 0, 0, 1]"^. When normalized, this eigenvector gives the state l^*) above. 



6 Construction of a state and observables which represent 
correlation functions 

In this section, we give a general construction, that is, states and operators which can 
represent all correlation functions. By using the theory of stabilizer states, especially, the 
binary representation of stabilizer formalism, we show that we can always choose a standard 
singlet for a state to represent all correlation functions. We believe this is the first time an 
explicit construction has been given. For the convenience of readers, we give a review of 
the binary representation of stabilizer formalism in [BJ 



We are interested in the construction of the states |^') and observables Ai and Bj 
satisfying \Ai\ < I, \Bj\ < I (or \Ai\ = \Bj\ = I), and 

mAi(S)Bj\'^) = {u'\v^) (14) 

in the case when two sets of vectors {u^}^^ and {v^Yj^^ satisfying the fourth condition 
(or the third condition) of Theorem [1] are given. As we have seen in the proof of the 

theorem, suppose ^ = min ^dimspanju*}™, dimspan{w-'}^^ for the fourth condition, and 
= d for the third condition. Then, these vectors are represented in a Hilbert space 7i 
whose dimension is 2*^ (i^ = | for even ^, and u = for odd ^) as Eq. (fT2|) by using 

a representation of the generators of the Chfford algebra and their simultaneous 

eigenvector \^). Therefore, in order to construct |^'), Ai and Bj, we need to know how we 
can find a simultaneous eigenvector l^') corresponding to In what follows, we give 

an explicit way to do it. 

First, we consider the case when ^ is even; that is, ^ = 2z^; In this case, we can choose 
the Weyl-Brauer matrices as an irreducible representation of the generators of the Clifford 
algebra on Ti.A and TCb SiS follows (see |A]) : 

fc— 1 u—k 

Xk 'Z (g) ■■■ (g) Z(g)X(g)'l(g)-'^--(g)l\ 

k—1 u—k 

Xk+u =^ -'Z (g -^^ ■ Z (gf^^^"^, (15) 

where X, Y and Z are the Pauli matrices, and an index k satisfies 1 < k < v [30j. Since 
{Xk}f!Li is a set of Hermitian operators which are composed of tensor products of Pauli 
operators satisfying Eq.©, {X^ ® Xk}'j!Li generate a commutable subgroup of the Pauli 
Group [15l[28l|T7] on 2v qubits 7^21^; we write this commutable subgroup as S2u- Note that 
Eg. ([15]) guarantees independence of the generator {XkfgXk^f^i- Since |^') is defined as the 
simultaneous eigenvector of {^fc®-^A:}fc=i corresponding to an eigenvalue 1 for all X^^Xk, 
\^) is nothing but a stabilizer state corresponding to the independent and commutative 
Pauli group elements {Xk (g X^}- Here, in a n-qubit Hilbert space, a stabilizer state |^') 
of commutative and independent Pauli group elements {M^.}^^^ is defined as a state which 
satisfies Mi \ ^) = \^) for all i |28l I17j . In this case, a set of 2" commuting Pauli operators 

S = {M £ Vn\M \^!) = 1^)} is called a stabilizer group of the state |^'); thus, {Mk}'^=i are 
independent generators of the commutative group S. Conversely, a set of 2" commuting 
Pauli operators all of which have real overall phase ±1 uniquely define a corresponding 
stabilizer state by the equations M\^) = \^) for all M £ S. Thus, our task is to find a 
concrete expression for a stabilizer state |^) corresponding to the stabilizer group S2u which 
is generated by {Xk g) Xk}l__-^. 

For the case when ^ = 2i/ + 1, we can choose an irreducible representation of the genera- 
tors of the Clifford algebra on TLa®'Hb as {^fc}^^^^ where {^fc}|=i are defined by Eq. p3]) 

u 

and X2U+1 '= Z (g ■ ■ ■ (g Z (see Since the equality X2U+1 = [—'iY^^=i^k holds, we 
derive the equality X2U+1 (gX2u+i = {—^Y^^=iXk ®Xk € S2u- Thus, a stabilizer state 1^*) 
of the stabilizer group S2U automatically satisfies X2y+i (g X2U+1 |^) = (— l)^+^fe=i |^). 
Therefore, by defining a2v+i = + Z]fc=i ^k, the state |^') is nothing but the state which 
satisfies (*| Ai ® Bj |^) = (n>^') for Ai and Bj defined by Eq.([l2]). 



So far, we showed that the state l^*) satisfying Eq. (|14p is a stabihzer state of the group 
which is generated by {X^ -^fc}fc=i- In order to construct a simultaneous eigenstate 
of {Xk ® Xk}1'Li defined by Eq. (fT5]) . we find a graph state which is equivalent to |^') under 
local Clifford operations by using the binary representation of the stabilizer formalism, 
which was mainly developed in [T5] and [25] (See[B]). Here, a graph state is a state in the 
standard form of stabilizer states [T7j, and defined by means of an n- vertex graph G with 
adjacency matrix {Oij}ij as the stabilizer state generated by commuting Pauli's operators 
{Kj}^^^ defined as 



Kj = x(^)n^=i ( z^'^) 



(16) 



Here X^'^ and Z^'^ are the Pauli operators having X and Z in the ith position of the tensor 
product and the identity elsewhere, respectively. We note that the binary generator matrix 
K corresponding to {iCj}"^]^ can be written as the following simple form: 



K 



def 



(17) 



In the binary representation, the set of Pauli operators {X^ ® defined in Eq.([T5]) 

is written as 

k-l v-k+l k-l v-k+l 
k—l u—k k—1 v—k 

1 I 1 o'^)^ 

k v—k k v—k 

Xk+u®Xk+^ ^ (1^0^ I 1^0^ I 

k—l v—k k—l v—k 

1 I 1 0^^)"^ , 

where \ < k < v. The generator matrix E corresponding to a set of Pauli operators 

\2v 

ifc=i 



{Xfc (gi Xfc}fc=i can be written as 



E 



def 



/ T,-i 


Tv 




Tu-i 






I 


I 




\ I 


I 


) 



(18) 



where T^^-i and Ty are v x u matrices defined as 

/ 1 1 
1 



def 



v-1 



def 





V 

/ 1 1 

1 



•• 1 \ 

• • 1 

1 
•• 0/ 

1 \ 

1 



ly 



and I is a v X 1^ identity matrix. 

In order to find a graph state which is equivalent to |^') under local Clifford operations, 
we try to transform the above generator matrix E into the form of Eq. (jl7p . which corre- 
sponds to a generator matrix of graph states, by changing the basis of the stabilizer and 
using local Clifford operations. As we can see in[Bl a change of the basis of a stabilizer cor- 
responds to applying an invertible matrix from the right-hand side of E, and an application 
of local Clifford operations to the state |^) is equivalent to applying a symplectic matrix 
from the left-hand side of E. Following a discussion which is similar to that used in [28j, 
we can find and easily check the following equation: 



L ■ E ■ R\ ■ R2 — F^i 



(19) 



where L, Ri, R2, F are defined as 



L 



def 



R'. 



dcf 



7 

/ 

0/0 
0/00 



/ 



F 



R 



def 



dcf 



/ 

/ 

/ 

/ 



Tu-1 






Tu 


I 


I 


I 


I 



(20) 



We also easily see that L is symplectic, Ri and R2 are invertible, and F is in the form 
of Eq. (jl7p : that is, F is a generator matrix corresponding to a graph state. Actually, in 
comparison with Eq. ijlTh . in this case, the adjacency matrix of the graph state corresponding 

^ ^ ' A way of constructing an arbitrary graph state is already known [Tj 



to F is 6* 



/ 



Thus, we can easily see the graph state corresponding to F is 







(21) 



where 1+) 



(|0) + |1)) and \-) = ^ (|0) - |1)) and that |*') is equivalent to n singlets. 
As we can see in[Bl in Eq. pop . applying Ri and R2 corresponds to changing the genera- 
tors of a stabilizer group of |^'), and applying L actually corresponds to applying a Clifford 
operation on 1^*). As we can easily check, a binary symplectic matrix L corresponds to a 

"local" Clifford operation £=^/(gi---(gi/(gi//(gi---(gi//, where the first n-qubits belong 
to TCa, tbe next n-qubits belong to TCb, and // is a Hardmer operator. Thus, the stabi- 
lizer groups corresponding to 1^*) and C"^ 1^*') = C |^') have the same binary representation 
E = ■ F ■ /?2^ • Ri^ ■ This fact does not guarantee l^*) = However, this fact 

guarantees that C |^'') is a simultaneous eigenvector of {Xi^(^Xi^}f,'Li, which is a stabilizer of 

l^*). In fact, we can easily check that £ |^'') = <S)i=i '^(\~^~^) ~^ I )) satisfies the following 

equations: For 1 < k < v, 



Xk+,0Xk+uC\^') = -C\^') 

X2,+l^X2,+lC\^') = C\^'), 



where we only use X-2u+i in the case ^ is odd. Therefore, as we have aheady mentioned, we 
can set = for \ < k < v and k = 2v + 1^ and = 1 for v + \ < k < 2^, and we also 

define Ai and Bj by Eq.l^^i). Then, j^f") =^ £ j^f') satisfies 

{^"\Ai(S)Bj\^") = {u'\v^). (22) 



Finally, as a result, by applying H ■ ■ ■ H to both |^'") and {X^ X^}'^'^^, we derive 
the following Theorem. 

Theorem 2 Given vectors {u"^}^^ and {v^}"^^ on Kinin(m,n) g^^igfyijig |^«| < |^i| < i 
for all i and j, there are operators {Ai}^^ and {Bj}'j^-^ on 2v qubits which satisfy \Ai\ < I, 
\Bj\ < I, and 

{^X''Ai^B,\^+f'' = {u'y), (23) 

where \^+) '= |00) + |11), = [^\ , C = niin ^dimspanju*}^, dimspanju-^}"^ , and Ai and 
Bj are defined as 

k=l 

Bj = Y.{-ir'{v'^)kXk. 

k=l 

In the above, n'* and v'^ are vectors derived by projecting and onto span{ti*}^-^ in the 
case o/dimspan{u*}™^ < dimspanjw-'}^^^ and onto spanju-'}^^ otherwise, 





def 


k~ 


1 


u—k 


Xk 






D I (g) • • • (g) I (1 < A: < 




dcf 


k- 


1 




Xk+u 






5 I® ■■■ (g) / (1 < A: < 1/ 




def 


V 






X21/+1 


X®-- 










= 


{l<k<iy 


and k = 2v + 1) 



and 



ak = I {u + l<k< 2v). 

That is, by choosing observables as above, n copies of a standard 2-dimensional maximally 
entangled state gives a quantum correlation which is represented by a correlation function 
{u^\v^). Although we used the fourth condition of Tsirelson's theorem (Theorem [1]) in the 
above theorem, we can also use the third condition of the theorem. In this case, we just 
need to modify ^ as ^ = d instead. 



7 Optimality of the construction 

As we can easily see from the discussion in the previous section, by choosing observables ap- 
propriately, the state 1$+)®*^ actually gives any quantum correlation function cij = {nJ'\v^') 



satisfying ^ < 2u + 1, where ^ = min ^dimspan{ti*}™, dimspan{u''}"j . Moreover, when 
a quantum correlation function Cij is an extremal point of the set of all m x n quantum 
correlation functions, this choice of a state and observables is actually the unique optimal 
choice with respect to both the dimension of the Hilbert space and the amount of entangle- 
ment, up to local unitary equivalence. This fact was also proved by Tsirelson ^27j. Before 
we state his theorem, we should note one important fact. For an extremal cij, ^ is uniquely 
determined, and called the rank of Cij [27\. That is, cij = (^u^\v^) = (^u'^\v'^) implies ^ = C' 
in this case. Here, we just give the theorem, when the rank ^ is even. The proof is found 
in [27!. 

Theorem 3 (Tsirelson) Let Cij be an extremal point, with even rank ^, of the set of all 
m X n quantum correlation functions. Suppose a state p and observables {Ai}'^^,{Bj}'j^^ 
on Ha Ti-B satisfy Eq.^. Then, Ha and Hb can be decomposed as Ha = ®a=i'^Aa 
and Hb = ®p=i '^Bfi such that dim'H^^ = dimT^^^ = 2^/^ for 1 < a < L and I < (3 < L. 
In this decomposition of the space, Ai and Bj can be written as 

Bj = (e^=isf )e/L+i, (24) 

and p can be written as 

L 

p = ^A,|$,)($„|, (25) 

a=l 

where {Xa}^^i is a probability distribution, and \^a) {^a\ is supported by Ha^ ^'Hb^- 
Moreover, the state \^a) and observables {A'^Y'j^^ and {B°'}^^ on Hac '^Ba local 
unitary equivalent to the state j^)*^^/^ and observables {Ai}'^^ and {Bj^^^^ defined in 
Theorem which corresponds to the same quantum correlation function Cij . 

Thus, the construction given in Theorem [2] is optimal with respect to the dimension of the 
Hilbert space. For odd ^, the corresponding theorem is slightly more complicated. However, 
the construction given in Theorem [2] is optimal also in this case |27J. Here, we note one 
more important fact which we can observe immediately from the above theorem. Since the 
observables {Aj}^^ and {Bj}'j^^ defined by Eq. (p^ obey the anti-commutation relation, 
they are actually proportional to a complex representation of the generators of the Clifford 
algebra C/(M^). Thus, a set of observables corresponding to an extremal quantum corre- 
lation function always has the structure of a Clifford algebra. Therefore, the construction 
using an irreducible representation of the Clifford algebra is optimal. 

We can easily calculate the value of entanglement measures for the state p defined 
by Eq. (j25p . (The reader unfamiliar with entanglement measures is referred to |23j.) 
First, since we can deterministically derive |$_|_)®^''^ by applying projective measurement to 
distinguish {HAcy}a=iJ derive a lower bound on the value of the distillable entanglement 
of Ed{p) > ^/2. Second, since Eq. (p5]) gives a pure-states-decomposition of p, we derive 
an upper bound on the entanglement formation Ep{p) < Moreover, since p defined by 
Eq. (j25p is maximally correlated, the entanglement formation Ef{p) is additive for p and it 
coincides with the entanglement cost Ec{p). Since any entanglement measure E{p) satisfies 
Ed{p) < E{p) < Ec{p), we can conclude E{p) = £;(|$)®^/^) = ^/2. Therefore, even if we 
use a mixed state p, we cannot reduce the amount of entanglement which is needed to 



represent an extremal quantum correlation function Cij, and we actually need ^/2 singlets 
to make such quantum correlation. In this sense, the construction given in Theorem [2] is 
an optimal construction with respect to both the dimension of the Hilbert space and the 
amount of entanglement. 

We conclude this section with the following remark. In some cases, the quantum state 
constructed above may be too large for any conceivable practical implementation. However, 
there is a way to approximate the given correlations using a state with much smaller di- 
mensions. This is by using the Johnson-Lindenstrauss Lemma, which we state in the form 
given in [9]: 

Lemma 1 (Johnson-Lindenstrauss) For e £ (0, 1) and N a positive integer, let K be a 
positive integer such that 

K > 4(eV2 - €^/?>y^log N. (26) 
Then for any set V of N points in there is a mapping f -.W^ ^ such that for all 

(1 - 6) {u\v) < {f{u)\f{v)) < (1 + e) {u\v) . (27) 

To use this lemma, we set = m+n, d = min{m, n), and let V be the union of the point 
sets {a^fclfcLi ^'^d {yi}f=i from Theorem [2j The lemma tells us that we may replace V by 
point sets in R^, which is a logarithmic reduction in dimension if m and n are comparable. 
Theorem [2] applied to these lower dimensional point sets gives a quantum setting where the 
number of qubits required is also reduced logarithmically, and for which the correlations 
can be made arbitrarily close to those desired. 

8 The maximum violation of quantum correlation vectors 

As we mentioned in Section[3l quantum correlation vectors are not characterized by Tsirelson's 
Theorem [H which gives necessary and sufficient conditions for quantum correlation func- 
tions. We saw that the theorem gives necessary conditions for quantum correlation vectors 
that can be combined with the no-signalling condition. These necessary conditions enable 
us to bound the maximum quantum violation of (general) Bell inequalities [2]. However, it is 
not clear to what extent the no-signalling conditions improve the bound given by Tsirelson's 
theorem. 

We denote the set of quantum correlation functions (resp. quantum correlation vec- 
tors) by MQ-g(m, n) (resp. Qcut(™'")) following the notation of [2]. The condition 3 
of Theorem [T] (±1 version) can be represented by the elliptope |10l section 28.4], which 
is well studied in combinatorial optimization. The elliptope e{G) of a graph G = (V, E) 
with n = |y| nodes is the convex body consisting of vectors z = (zij) S such that 
there exist unit vectors u^,u^ S SR" for each node i,j G V satisfying Zij = (n*|n-^). Let 
Km,n be the complete bipartite graph with nodes Vm,n = {^i, •••! Am,Bi, Bn} and edges 
Em,n = {^iBj I 1 < i < m, 1 < J < n}. The elliptope e{Km,n) is the set of the vectors in 
[—1,-1-1]™" satisfying condition 3 of Theorem [H that is: 

Mq^{m,n)=e{K^^n)- (28) 
We can apply condition 3 of Theorem [T] to quantum correlation vectors as follows: 



• For a quantum correlation vector x = {cki) G there exist three sets 

of unit vectors {t^}h=i, {^^lilLu {^'■'}?=i £ I^^j and d < 1 + m + n, for which 



X = (cm) 



= i{t'\u')l<i<mAt'\v')i<,<. 
(^*l'^"')l<j<m,l<j<n)- 



The set of the vectors satisfying the latter conditions corresponds to the elliptope of another 
graph VKm,n called the suspension graph of Km,n- The graph VKm,n = (Vy^.n) V£^m,ri) 
is defined as the graph with 1 + m + n nodes Wm,n = {X} U Vm,n and m + n + mn edges 
'VEjn,n = {X^i I 1 < ^ < ""T-} U {^^j | 1 < J < U Em,n- The no-signalling conditions 
correspond to the rooted semimetric polytope denoted by RMet(VKm,n), which has been 
also well studied [101 section 27.2], and is given by the following inequalities in the ±1 
version: 



It has recently been shown that the inclusion in equation (j29p is in fact proper. Navascues, 
et al. [21] recently succeeded in completely characterizing the set of all quantum correlation 
vectors by means of a sequence of SDPs [21]. Using this characterization they got a tight 
upper bound on the quantum violation of the Froissart Bell inequality, 13322 (see [13], [8]) of 
.25089. This is significantly tighter than the bound of 0.3660 given by the right hand side 
of (I29p for this inequality [2j. A similar computational result was reported by Doherty et 
al. [11], also using a hierarchy of SDPs. The lower bound on the quantum violation is .25 



We begin by reviewing the work of Navascues, et al. First, they treat a more general 
problem setting. That is, their result is valid for all quantum correlation experiments which 
consist of measurements having any number of outcomes. So, we need to generalize the 
definition of a quantum correlation vector. Suppose Alice's and Bob's measurements are 
described by projective measurements {Ea^}"^^^ and {-^^"'^l^^p where a{i) and b{j) are the 
number of outcomes of the zth and jth measurements of Alice and Bob, respectively. That 
is, this set of operators satisfies 



XXA,+XxB,+XA,Bj > -1, 

XX Ai + XX Bj - XAiBj < 1, 

XXAi - XX Bj + XAiBj < 1, 

XXA,+XxBj+XA,Bj < 1- 



Quantum correlation vectors are bounded as follows [21 Theorem 6]: 



2Cut("^'^) ^ e{^K„,,n) nRMet{VKm,n)- 



(29) 



a(i) b(j) 




0=1 /3=1 



[i^«,4^')] = 0. 




(30) 



Then, by defining P(a|i) = {^\ E^a^ \^), P{P\j) = {^\ Ef \^), P{a,P\i,j) = (*| E^^^ E^f |^), 
a quantum correlation experiment is completely described by the vector 

{{P{a\i)}a,^ , {P{m}pd , {P{a, P\h j)U|3,^,j ) , (31) 

for 1 < a < a{i) — 1, 1 < i < m, 1 < f3 < b{j) — 1, 1 < J < n. The set of these vectors 
is denoted as Q{m,n,a{i),b{j)). We can easily see that Q{m,n, a{i),b{j)) is isomorphic to 
QQ^^(m, n) when a{i) = b{j) = 2 for all i and j. Suppose T is the set of all projections /, 

Ea^ and F^"'^ for all (a, /3, except {a{i), b{j), that is, a runs from 1 to a{i) — 1 and /? 
runs from 1 to b{j) — 1 . We represent a product of projections in T as a sequence on T (e.g. 
Ea\ Ea'^ E^,\ Ea'^ E^^i'^ F^"* are sequences on T). For a sequence S, we define the length 
of S as the minimum number of projections by which we can write S] as a convention, we 
define the length of / as 0. Then, we define Sc as the set of all sequences whose length is 
no greater than c; of course, we have Sc ^ Sc+i- Let Sc = {•Sfclit^i) where \Sc\ is a size of a 
set Sc- Then, we define a finite set of independent equalities J-{Sc) as a set of equalities of 
the form 

J2{F)M{'^\SkSi\^)=g{P) (32) 

M 

which is satisfied by sequences in Sc for an arbitrary choice of projections {Ea^}"^^^ and 
^pU) }6W satisfying Eq. §0^. InEq.^Si), g (P) is an affine function of probabilities P{a, f3\i, j) 
defined as 

g{P) = {g)o + Yl i9)a,|3,^,JP{a,^3\i,j). (33) 

J^{Sc) is not uniquely determined from this definition. But, for all c , there actually exists 
an easy algorithm to find a set J^{Sc) satisfying the above condition [20] . Since each equality 
in J-{Sc) is determined by a \Sc\ x \Sc\ matrix F and an affine function k on rI'^'^I, we write 
(F, k) e J^(5c) if F and k satisfies Eq.([32]) in J^(5c). 

Now, we are ready to present the main result of Navascues et al. [21]. Suppose 
Qc{m, n, a{i), b{j)) is defined as a set of quantum correlation vectors satisfying the following 
condition: there exists a \Sc\ x \Sc\ positive semidefinite matrix T > satisfying 

Tr {F^T) = g (P) (34) 

for all {F,k) G J^{Sc). We have Qc+i{m.,n,a{i).,b{j)) C Qc{m,n, a{i),b{j)) Their main 
result is that 

Q{m,n, a{i),b{j)) is actually a limit of the sequence of the sets {Qc("i, n, a(z), 6(j))}^Q, 
that is, 

oo 

Q{m,n,a{i),b{j)) = Q Qc{m,n,a{i),b{j)). (35) 

c=0 

Moreover, by the definition, a quantum correlation vector P is in Qdm, n, a{i), b{j)), if and 
only if 

max{A \T - XI > 0, 

Tr (F^r) = g (P) , V(F, k) G J^{Sc)} > 0. 

Since the above optimization problem is an SDP, we can check whether P is in Qdin, n, a{i), b{j)) 
by an SDP of the size \Sc\ x \Sc\- Thus, we can judge whether P is in Qc{m, n, a{i), b{j)) by 



this infinite sequence of SDPs. In [2T| , they also gave a condition such that if it is satisfied by 
r corresponding to a positive A, then we can immediately conclude P £ Q{m, n, a{i), b{j)), 
without having to solve an infinite sequence of SPDs. 

Navascues et al.'s method can be applied to the calculation of maximal violation of 
Bell inequalities. That is, we can always derive an upper bound of the maximal violation 
by optimizing Bell inequalities over Qc{m,n,a{i),b{j)), and this optimization is also SDP 
[2T]. Suppose we are interested in a Bell inequality J{P) = Y2af3ij ^(Q^>/55^)i)-P(o!,/3|i,j), 
for some coefficients V{a, j). Then, the optimization of J{P) over Qc{m,n,a{i),b{j)) 
reduces to the following SDP: 

max{J(P)|P E Q,{m,n,a{i),b{j))} 
= max{Tr(Ccr) | T > 0, 

and Tr {F^T) = g (P) , V(F, k) G T{Sc)}. (36) 

In the above equation, the \Sc\ x \Sc\ matrix d^c is defined as follows: Since each index of 
the matrix is associated to a sequence in Sc, we directly represent each index by the 
corresponding sequence. Then, (Cc)p(i) ^u) = ^ict, (3,i, j)/2 and {Cc)ss' ~ ^ other 

combination of sequences {S,S') G 5c x Sc- By Eq. ([35]) . we can easily see that Eq. ([36]) 
actually converge to the maximal violation of J{P)- 

Also very recently, Doherty et al. derived a sequence of SDPs which converge to the 
maximal violation of a Bell inequality pJLj. In this final part of this section, we review their 
result. Suppose P is a collection of Hermitian polynomials defined as follows: As variables 
we consider {Ea^}a,i and {Fp^}^^. The number of variables is equal to the number of 
projections which appear in the quantum correlation experiment, and are labeled in the same 
way as the projections. We should note that {Ea^}a,i and {Pg"'^}/?^ are non-commutative 
variables, and do not necessary satisfy condition Eq. (j30p . We define a set of polynomials 
Ti, T2, 73, T4, 75 on non-commutative variables {Ea^}a,i and {Fp^}i3j as follows: 

1 < a < a(i), 



1 < /? < b{j), 1 < i < m, 1 < j < n| 



^2 = (Uu-E^i^^}lu|Uu-E^^^} 



rs = \[j{{E, 



(^)^2 




U {f<''F5;> + f<?'F5'>}^^^,.. (37) 

P»1 



We define T = Ti U T2 U T3 U 7^ U T5, and P = T U -T. Then, for a Hermitian {E^a'ja i and 

/3 



{fI"'^}^ j, P is actually a set of Hermitian polynomials. We should note that values of all 



polynomial in V are zero for and {Pfl^^}/3j satisfying Eq.(l30]). Then, we define the 



convex cone C-p^c as follows: q is in Cp^c, then g is a polynomial on {Ea }a,i and {F^ }i3j 
whose order is no greater than 2c, and which has the form q = rjrj + J2i ^IjPi^ij^ 
where pi £ V, and r^, Sij are arbitrary polynomials on {Ea^a,i and {pj^-'^pj. Suppose a 
polynomial on {Ea^}a,i and is defined as 

q, = i^I-l ^(«' /3> hj)E^a^4'^ I • (38) 
Then, the optimization problems 

Wc =^ min{z^|gi, e Cp,c} (39) 

converge to the maximal violation of the Bell inequality J(P) in the limit c ^ oo. Moreover, 
surprisingly, Wc can also be calculated by SDP [llj. So, this optimization also gives another 
way to calculate an upper bound of maximal violation of a Bell inequality as close as we 
like by SDP until the limit of our computational power. 

Finally, we note one important fact. Even though Q{m,n,a{i),b{j)) can be character- 
ized by a sequence of SDPs, this does not imply that Q{m, n, a{i), b{j)) can be characterized 
by a polynomial time algorithm. For example, if the convergence of Qc{m,n,a{i),b{j)) to 
Q{m,n, a{i),b{j)) is slow, then Eg. ([35]) does not guarantee the existence of such a polyno- 
mial algorithm. Actually, we can see a similar situation in the characterization of entangled 
states. The set of all separable (non-entangled) states can be characterized by a sequence of 
SDPs. However, the problem to check whether a given state is separable or not is actually 
NP-hard It is an important and interesting open question to see whether or not a 

similar result holds for quantum correlation vectors. 
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A Clifford algebra 

In this appendix, we give several important properties of the Clifford algebra associated to 
a real linear space V = M", which are used in this paper. All of the facts given in this 
appendix can be found in p^l [30]. 

The Clifford algebra Cl{V) associated to a real linear space V = M" is defined to be the 

quotient algebra T{V)/I{V), where T{V) =^ is the tensor of V, and I{V) is the 

ideal of T{V) generated by all element of the form v (8) v— < v\v > 1 for v £ V. Although 
this is the definition of Clifford algebra, the following well-known characterization is more 
convenient in this paper. 



Theorem 4 Let ei, • • • , e„ be any orthonormal basis of a real linear space V , then, Cl{V) 
is generated 6y ei , • • • , e„, subject to the relations 

CiCj + CjCj = 25ij. (40) 

Thus, a set of aU elements in the form Cj^ • • • with < • • • < i/j is a basis of Clifford 
algebra Cl{V), and the dimension of Cl{V) is 2". 

In this paper, we are interested in the ^^complex representation" of Clifford algebra 
Cl{V), which is a real algebra homomorphism from CliV) onto an algebra of all linear 
operators ^{7i) on a complex finite-dimensional Hilbert space 7i. The following property 
of irreducible complex representations of CKM"") is important for our paper. 

Theorem 5 For even n, up to equivalence of representations, Cl{W^) has only one com- 

n 

plex irreducible representation, whose dimension is 2 2", and which is faithful. For odd n, up 
to equivalence of representations, CZ(M") has only two different complex irreducible repre- 

n-l 

sentations, whose dimensions are 2 2 , and which are not faithful. 

We can make these irreducible representations in the following way. Suppose n = 2j^ is 
even, and W = is a z^-qubit Hilbert space. Then, on Ti., the Weyl-Brauer matrices are 
defined as 

i—1 v—i 

i—1 u—i 

Xi+^ =^ -'Z ■ Z (^Y ^7^^^"^, (41) 

where X, Y and Z are the Pauli matrices, and an index i satisfies \ < i < v. We can 
easily check that these matrices satisfy the anti-commutation relation (j40p . Thus, they 
are a representation of the generators of Cl{W^). Since they also generate the algebra of 
all linear operators !B('H) which has the same dimension of C/(M"), this representation is 
faithful and irreducible. For odd n = 2u + 1, we can choose an irreducible representation 
of the generators of the Clifford algebra Cl{W^) on a n-qubit Hilbert space 7i as {^iji^i"^ 

V 

def / ^ \ 

where {Xi]'^^^ are defined by Eq. (fT5]) and X2v^i = Z ® ■ ■ ■ ® Z . Again, we can check 
that the anti-commutation relation (j40p . and {Xi\^^^^ generates the algebra of all linear 
operators !B('H), which agrees with the irreducibility of this representation. However, this 
representation is not faithful, since now the dimension of this algebra is smaller than that 
of CZ(M"). Another inequivalent irreducible representation is derived by redefining —Xi as 
Xi for 1 < i < u. 

B Review of the binary representation of the stabiUzer for- 
maUsm 

In Section [6] of this paper, we use the binary representation of the stabilizer formalism. For 
readers who are not familiar with this topic, we give a small review of it. All the detailed 
proofs of the statements that appear in the following part can be found in [15] and |28j . 

First, we rewrite the Pauli matrices by using binary indices as / = aoo,-'^ = foi,^ = 
cjii,Z = cJiQ. Then, this notation gives an encoding of the Pauli matrices into the binary 
field Z2 as 

I ^ 00, X ^ 01, Y ^ 11, Z ^ 10. (42) 



More generally, we can encode the tensor products of Pauli matrices into a vector space on 
the binary field as: 



def 

Suppose P is a 2n X 2n binary matrix defined as P = 



where / is an n x n 



def fc>, M 

^ {ui ■ ■ ■ un\vi ■ ■ ■ Vn)'^ , (43) 

where u =^ (ui ■ ■ ■ un), v *== (vi - ■ ■ V]\f), and T is transposition. We can easily see this 
representation satisfies o"(^^^)(T^^, ^-j^ = aa^^^^, ^^^y where a is a phase factor which ranges 

over {±1, ±1} depending on u, v, u',v' , and © is bitwise summation in the binary field. Since 
the Pauli group on n-qubits Vn can be written as 

Vn = {a(T^u,v)\a G {±l,±i},u,v eZ^}, (44) 

Eq. (j43p gives a representation of the n-qubits Pauli group on the 2n dimensional vector 
space on binary fields Z^"; this representation is not faithful since we omit the phase factor 
a in Eq. ([H]) . We should note that the correspondence ([13|) can be also understood as giving 
the projective representation of the finite field on n-qubits Hilbert space. 

/ 
/ 

identity matrix. Then, "(T(^^^) commutes with a^^, if and only if {u,v)P{u' ,v')^ = 0." 
Thus, commutativity of elements of the Pauli group is equivalent to the orthogonality of 
the corresponding binary vectors under the symplectic inner product; {u,v)P{u' ,v')'^ is 
called the symplectic inner product between {u,v) and {u',v'). Further, suppose 5 is a 
commutative subgroup of the Pauli group on an qubit system. Then, we can easily see 
the following fact: {Xi}^^ are m independent generators of S if and only if a set of vectors 
{^HfcLi on corresponding to {Xjj^j^ is an orthogonal basis of a subspace S of 
corresponding to S in terms of the symplectic inner product; in other words, the generator 

matrix E '= (el, • • • ,6^) satisfies E^PE = 0. In this case, dimS = m and all pairs of 
different vectors e and / on © satisfy e^Pf = 0. In particular, the stabilizer group S of an 
A^-qubit stabilizer state corresponds to an dimensional self-dual subspace S; that is, S 
is equal to its symplectic orthogonal complement S"^ = {e G Z^^le^P/ = (V/ S S)}. 

Although a stabilizer group of a stabilizer state determines a corresponding self-dual 
subspace, a given A^-dimensional self-dual subspace of a binary field, or a given generator 
matrix does not uniquely determine a corresponding stabilizer group and a stabilizer state. 
This fact can be seen as follows: Suppose {Mj}^^ is a set of commutative and independent 
Pauli operators on A/^-qubit Hilbert space, and E \s a. generator matrix corresponding to 
{Mi\f^^. Then, for all (ci, •••,€„) G Zjv, {{—'^Y^ Mi}f^^ has the same matrix E as its 
corresponding generator matrix. In this case, although a stabilizer state |^') corresponding 
to {{—ly^ Mi}f^-^^ is not a stabilizer state corresponding to {Mi\f^^, \^) is a simultaneous 
eigenvector of {Mi]f^^. 

Suppose {ek}^^i is an orthogonal basis of a self dual subspace Sat, and suppose R is 
an A^ X A^ invertible binary matrix. Then, a set of vectors {fi}iLi defined by 

N 

r def 



fi = Y.e-lRu (45) 



k=l 



is another orthogonal basis of Sat. Suppose {X^}^^^ are independent generators of a 
stabilizer group Sm, {el} are their binary representation, and P is an A^ x A^ invertible 



binary matrix. Then, {afc^A:}fc=i is another set of independent generators of Sn, where 
Ofc G {ibljibi} is a suitably chosen phase factor, and is an element of the Pauli group 
corresponding to a vector fk defined by Eq. (jl5]) . Note that if corresponds to the binary 
vector /fc, then, ±1^, and izYfe also correspond to fk- 

The Clifford group ^at is defined as the set of all unitary operators satisfying UVnW C 
Vn foi" the qubits Pauli group Vn- We can easily see that an element U of Clifford 
group induces the linear transformation fu : — > Z^^. In this case, it is known 
that the set of such transformations {fu}uegN coincides with the set of all symplectic 
transformation on Zg^. Thus, for an element of the Clifford group U, Ua^U^ = aa^, for a 
suitably chosen a £ {±1, ±i}, if and only if there exists a symplectic matrix Qjj (a matrix 
satisfying QjjPQu = P) such that /[/(e) = Que = e'. 

C Computational investigation of the strength of the no- 
signalHng inequahties 

We may wonder to what extent the inclusion (j29p is essentially stronger than Tsirelson's 
construction. In this appendix we review an investigation of a set of known Bell inequal- 
ities computationally, showing that in most cases a tighter bound on quantum violation 
is achieved by using the no-signalling conditions. In some cases, however, the bounds are 
identical. 

Takahashi et al. [25] report on a computational investigation of the optimal value of 
Bell inequalities with the following constraints: 

(i) the elliptope e{VKm,n) 

(ii) the elliptope eiVKm.n) and the rooted semimetric polytope RMet(Vi^m,n) 

The optimal values were calculated by the method of Avis, Imai and Ito [2] using the software 
SDPA [14] for a set of 89 Bell inequalities [19], which includes two quantum correlation 
functions: A2 (the CHSH inequality) and A8 (the 13322 inequality [8]), and the other 87 
quantum correlation vectors. The results are reproduced here as Table 1. As A2 and A8 
are quantum correlation functions, it is consistent with Tsirelson's Theorem [1] that both 
optimal values are the same. On the other hand, we also see that for two Bell inequalities, 
A19 and A56, both optimal values coincide. This means that, for the Bell inequalities A19 
and A56, the rooted semimetric polytope RMet(V^Cm,n) is inactive. 
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Table 1: The values of optimum solutions of the semidefinite programming problems for 
those 89 Bell inequalities [19j 
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